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Abstract. We prove the local in time existence of the classical solutions to the 
system of equations of isothermal viscoclasticity with clamped boundary conditions. 
jyy \ We deal with a general form of viscous stress tensor Z(F, F), assuming a Korn-type 

condition on its derivative DpZ(F, F). This condition is compatible with the balance 
of angular momentum, frame invariance and the Claussius-Duhem inequality. We give 
examples of linear and nonlinear (in F ) tensors Z satisfying these required conditions. 



1. Introduction and the main results 

In this paper, we are concerned with the local in time existence of the classical 
solutions to the system of equations of isothermal viscoelasticity. The system we study 
is given through the balance of linear momentum: 

(1.1) &-div(w(V0 + 2(V£,V&)) =0 infixM+! 
and it is subject to initial data: 

(1.2) £(0,-) = £o and 6(0, •) = 6 in fi, 
the clamped boundary conditions: 

(1.3) £(-,x) = x vx e an, 

and the non-interpenetration ansatz: 

(1.4) detV£>0 in a 

Here, £ : Q x M + — y IR n denotes the deformation of a reference configuration Q C R™ 
which models a viscoelastic body with constant temperature and density. A typical 
point in Q is denoted by X, and the deformation gradient, the velocity and velocity 
gradient are given as: 

F = V£ G R nxn , v = & e K n , Q = V& = Vv = F t e R nxn . 

In (11. ip the operator div stands for the spacial divergence of an appropriate field. We 
use the convention that the divergence of a matrix field is taken row-wise. In what 
follows, we shall also use the matrix norm \F\ — (ti(F T F)) 1 ^ 2 , which is induced by the 
inner product: F\ : F 2 = tr(F 1 T F 2 ). To avoid notational confusion, we will often write 
(Ft : F 2 ) instead of F ± : F 2 . 
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1.1. The elastic energy density W. The mapping DW : IR nxn — ► M nxn in ( 11. ip is 
the Piola-Kirchhoff stress tensor which, in agreement with the second law of thermody- 
namics [8], is expressed as the derivative of an elastic energy density W : M nxn — > M + . 

The principles of material frame invariance, material consistency, and normalisation 
impose the following conditions on W, valid for all F G M nxn and all proper rotations 
R G SO(n): 

(i) W(RF) = W(F), 

(1.5) (ii) W{F) -> +oo asdetF^O, 
(hi) W(Id) = 0. 

Examples of W satisfying the above conditions are: 

W\{F) = \{F T F) 1/2 -Id| 2 + |logdetF| 9 , 

W 2 (F) = |(F T F) 1/2 -Id| 2 + 

where q > 1 and is intended to be +oo if detF < [22]. Another case-study 
example, satisfying (i) and (iii) is: W (F) = \F T F - Id| 2 = (\F T F\ 2 - 2\F\ 2 + n). 

We will assume that W is smooth in a neighborhood of SO(n). Since div(LW(V£)) 
is a lower order term in ( 11. ip . it follows that other properties of play actually no 
role in the proof of our main Theorem 11.11 and 11.21 We hence remark that the same 
results are valid when dw(DW(V^)) is replaced by div(Diy((V^)A(X)^ 1 )). Such term 
corresponds to the so-called non-Euclidean elasticity, where the deformation £ of the 
reference configuration strives to achieve a prescribed Riemannian metric g = A T A on 
Q. This model pertains to the description of prestrained materials and morphogenesis 
of growing tissues [191 HH] • 

1.2. The viscous stress tensor Z. The viscous stress tensor is given by the mapping 
Z : W nxn x R raxrt — y ]R nxn , depending on the deformation gradient F and the velocity 
gradient Q. It should be compatible with the following principles of continuum me- 
chanics: balance of angular momentum, frame invariance, and the Claussius-Duhem 
inequality [8]. That is, for every F,Q G IR 3x3 with det F > 0, we require that: 

(i) skew (F _1 Z(F, Q)) = 0, i.e. Z = FS with S symmetric. 

(ii) Z(RF, R t F + RQ) = RZ(F, Q) for every path of rotations R : R + — ► 

(1.6) SO(n), i.e. in view of (i): S(RF, RKF + RQ) = S(F, Q) WR G SO(n) 
\/K G skew. 

(iii) Z(F,Q) : Q > 0, i.e. in view of (i): S : sym(F T Q) > 0. 
Examples of Z satisfying the above are: 

Z m (F,Q) = [sjmiQF-^+'F-^, 

(1.7) Z' (F, Q) = 2(detF)sym(QF- 1 )F- 1 < T , 

Z^F,Q) = 2Fsjm(F T Q). 



detF 



for det F > 0, 
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We note that in the case of Z' Q , the related Cauchy stress tensor Tq = 2(detF)~ 1 Z2F T = 
2sym(QF~ 1 ) is the Lagrangean version of the stress tensor 2symVw written in the 
Eulerian coordinates. For incompressible fluids 2div(symVt>) = Av, giving the usual 
parabolic viscous regularization of the fluid dynamics evolutionary system. 

1.3. The main results. Our main assumption implying the dissipative properties of 
(11. ip will be expressed in terms of the following condition on a (constant coefficient) 
linear operator M : R nxn M nxn : 

(1-8) \\V(\\l 2m < 7 f (MVC) : VC VC G W%{R n ,R n ). 

Note that (11.81) is a Korn-type estimate, reducing to the classical Korn inequality for 
M(F) = symF and 7 = 2 [16]. Naturally, (II. 8p is equivalent to (12. ip which is the same 
estimate valid for all C e W?(U, R n ) with (\au = 0, on any fixed open, bounded U C R n . 
It can be shown, via Fourier transform (see Lemma |2~2]) , that (II. 8p is also equivalent to 
the strict positive definiteness of M when restricted to the space of rank-one matrices 
Q = a ® b: 

(1.9) Va,beR n \a\ 2 \b\ 2 = \a®b\ 2 < 7 (M(a <g> 6) : a® 6). 

We point out that the above condition resembles, naturally, the local well-posedness 
criterion for the inviscid elasticity system [15], where the validity of (11.91) for M = 
DW(V£,o) is equivalent to the hyperbolicity of the first-order system (II. ip with Z = 0. 

The main result of this paper is the following: 
Theorem 1.1. Let Q be a smooth bounded domain in M n and let: 

(1.10) £0 g W 2 (n), 6 e W 2 ' 2/p (n) for some p>n + 2, 
satisfy: 

inf det V£ W > 0, £ (X) = X, ^(X) = Vie dQ. 

Assume that the viscous tensor Z has the property that: 
(1.11) 

\/X G n, / TOP holds with M = DqZ(V£ (X), V^(X)) and 7 independent of X. 

Then there exists T max > such that the problem M.l\) , U.ty) , M.3\) , fli.^| ) admits a 
unique regular solution £ G W 2 > 2 (Q x (0, T)) fl Loo((0, T), W 2 (Vt)) with & G W 2 ' l (Q x 
(0,T)) for all T < T max . 

The proof of Theorem ll.il will be given in sections [2J [31 HI In section Owe show that 
viscous stress tensors in (I1.7P satisfy (II. lip : for any initial data £0, £1 in case of the linear 
in Q tensors Z , Z' 0) Zq, and for initial data enjoying additionally det sym(V^i(V^o)~ 1 ) 
in case of the nonlinear (in Q) tensors Z\,Z2, see Lemma 12.31 Thanks to this ob- 
servation, Theorem 11.11 proves the mathematical well-posedness of a class of physically 
well-posed models. 
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With the same techniques of proof of Theorem 11.11 one can show that: 

Theorem 1.2. Let S be the solution operator of the problem U.l\) - ljl-4\ ) described 
in Theorem \l.l\ given by: 

5:^ 2 (0)x^(n) (^- 2 (fi x (0, T)) n ^((0, T), lf p 2 (fi))) x x (0, T)). 

Then S is continuous. 

We omit the proof and refer instead to standard texts [TJ [171 EH EI], or to an 
application of the same methods in the more current context as in Theorem 1.2 [9] . 

1.4. Relation to previous works. The dynamical viscoelasticity (11. ip has been the 
subject of vast studies in the last decades. For Z(F, Q) = Q conflicting with the frame 
invariance (jl.6p (ii), various results on existence, asymptotics and stability have been 
obtained in (2j EH EH [12]. For dimension n — 1, existence of solutions to ( 11. ip has 
been shown in [7J H] for Z depending nonlinear ly on Q. 

Existence and stability of viscoelastic shock profiles for a large class of models orig- 
inating from (II. ip has been studied, among others, in [31 [5] . 

Existence of Young measure solutions to system ( II. ip was shown in [TO], without any 
additional assumptions on Z, but with condition ( II. 6p (iii) strengthened to the uniform 
dissipativity i.e: Z(F,Q) > 7|Q| 2 . These measure-valued solutions were shown to be 
the unique classical weak solutions under the extra monotonicity assumption: 

(1.12) (Z(F 1 , Q 1 ) - Z(F 2 , Q 2 ) : Q 1 - Q 2 ) > K \Q X - Q 2 \ 2 - l\F, - F 2 \ 2 , 

see also [25] for a treatment of slightly more general type of PDEs under the same 
condition. As noted in [10], (I1.12p is incompatible with the balance of angular mo- 
mentum (ll.6p (i). In particular, (I1.12p is not satisfied by any of the examples in ( II. 7p . 
even Z , Z' Ql Zq which enjoy condition (jl.lip for any invertible F = V£oPO an d an Y 
Q = VCi(X). 

From the theory of PDEs viewpoint, our present result is a rather straightforward 
application of the theory of nonlinear (quasilinear) parabolic systems. Namely, we 
apply the maximal regularity estimates to control the nonlinearities of the system (11.11) . 
We choose the L p -framework in order to avoid technical difficulties, but a similar results 
and estimates are expected in the Besov spaces framework [9]. In a sense, our result is 
hence a consequence of the classical works of Ladyzhenskaya, Solonnikov and Uralceva 
[17] , which has been further developed in [TJ [TTJ [20], and which is a powerful tool in 
the study of the parabolic-elliptic systems. 

1.5. Notation. By L P (Q) we denote the space of functions integrable with respect to 
the Lebesgue measure, with p-th power. By x (0,T)) for k, I G N we denote 
the anisotropic Sobolev space defined by the norm : 

W u \\w£> l (nx(o,T)) = IMk(nx( 0i T)) + l|VV<^w|| Lp( nx ( o,T)), 
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where V fc is the k-th space derivative and d t is the time derivative. The isotropic 
version is given by: 

The space Wp P (Q) is the trace (in time) space of Wp' x {Q, x (0, T)). For further 
details we refer to j6]. 
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2. The constant coefficient problem 

The following auxiliary result will be needed in the proof of Theorem 11.11 

Lemma 2.1. Assume that M : M nxri — > W axn is a linear map satisfying the Korn-type 
inequality: 

(2.1) ||VCHL(a) < 7 / (MVC) : VC VC G W%{U,R n ) with C\ 9U = 

Ju 

Then the solution to: 

( Ct-div(MVC) = / »f/x(0,T), 

(2.2) < C = ondU x (0,T), 

[C(0,-) = Co inU 

admits the following maximal regularity estimate: 

(2-3) \\Ct, V 2 C||l p (C/x(0,T)) < C^uiWfWLpiUx^T)) + HColl^a-a/p^), 

where the dependence of C on U is uniform for any family of domains which are 
uniformly bilipschitz homeomorphic to each other after appropriate dilations. 

Towards a proof of Lemma 12.11 note first that for M = Id, i.e. when (12. ip holds 
trivially, (I2.3P is a classical maximal regularity parabolic estimate for the heat equation. 
When M(F) = symF, i.e. when (12. ip reduces to Korn's inequality, the proof of (12. 3p is 
also immediate. For, take div of the equation in (12. 2p . and note that div T div(symVC) = 
div (|AC + |VdivC) = |(divAC + |AdivC) = AdivC so that: 

(divC)t - A(divC) = div/ in U x (0, T). 

By the maximal regularity estimate for the heat equation: 



(2.4) ||VdivC|U P (t/x(o,T)) < C pJU (\\f\\ LA ux(o,T)) + lldivCol 



w}- 1/p {u) 
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Now, (12.21) can be written as: 



Ct-^AC = /+^VdivC mUx (0,T). 

We hence obtain: 

IK*, V 2 C||l p (c/x(o,t)) < C P)U (\\f + ^VdivC|| Lp (c/x(o,T)) + ||divCo|| w 2-i/ P(c/) 

which combined with (12.4ft yields ( 12.31) . 

In the general case, Lemma |2~T1 follows from the maximal regularity theory developed 
for parabolic initial-boundary value problems in [IT]. Under the ellipticity condition 
(b) on page 98 in there (see also Definition 5.1), the estimate (12. 3p is a consequence of 
Theorem 7.11. We now prove that condition (12. ip implies that the constant coefficient 
operator — div (MVC) has its spectrum contained in the proper sector of the complex 
plane, which immediately gives ellipticity in the sense of [TTj . 

Lemma 2.2. Conditions U.8\) . \2. 1\) and U.9\) are equivalent. Moreover, under any 

of these conditions the operator —div MV(-) is elliptic, i.e: 

(2.5) 

spec(— div MV(-)) C u G C : Re z > 0, and arg z G [a*, a*] with < a* < a* < — >. 

L 2 2 J 

Proof. 1. Conditions (11. 8p and (12. ip are equivalent in view of the density of C^°(W n ) 
in Wi(R n ). To include (II. 9p . we use linearity of Fourier transform and Plancherel's 
identity: 



I|vcI|l 2( ^) 2 = II(vC) a IU 2 ^) 2 = / \((k)®k\ 2 dk 



(M(VC) : VC) = / (M((VC) A ) : (VC) A ) = / (M(((k) ® k) : (((k) ® fc)> d*. 
Hence, (jl.8p is equivalent to: 

(2.6) VC E W?(R n ) ||C ® fc||i 2(R „) < 7 / (M(C ® fc) : (C ® k) ) dk. 

It is therefore clear that ( II. 9p implies (ll.8p . On the other hand, given ko, ao G M n , con- 
sider: ( m (k) = (^pU, 2 (k — k ) + pU 2 (k + fc )J ao) where p m is the standard radially sym- 
metric mollifier supported in the ball B(0, 1/m). Applying (12. 6p to £ n e W 1 2 (R n , M n ) 
and passing to the limit m — > oo, yield (jl.9p for the matrix Q = oq ® fco. 

2. To prove (12. 5p . consider the eigenvalue problem: 

AC - div(M(VC)) = in M n , 
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which after passing to the Fourier variable k el B becomes: 
(2.7) \((k) = M(((k) <g> k)k. 

Upon writing A = cr\k\ 2 , the problem (12 .7p is equivalent to locating the eigenvalues a of 
the family of linear operators {M fc }| fc | =1 , M k : W 1 — > R n given by: M k (a) = M(a (g) k)k. 
Recalling ( II. 9p we see that each M k is strictly positive definite: 

\k\ 2 

M k (a) ■ a = (M(a <g> k) : (a g> fc)) > — lal 2 . 

Consequently, spectrum of every satisfies Re ex > 0. By continuity with respect to 
k which varies in the compact set \k\ = 1, we obtain the inclusion (12. 5p . ■ 



Finally, we have the following: 

Lemma 2.3. TTie viscous stress tensors Z in (fJ.Tp satisfy A2.1\) with M = DqZ(F , Qq), 
for every F , Q with det F > 0, in the following manner: 

(i) Zg with 7 = IFq" 1 ^! 2 . 

(ii) Z' Q with 7 = (Fol^detFo)" 1 . 

(iii) Z G with 7 = \\F \ 2 . 

If we additionally assume that det sym(Qo-^o 1 ) 7^ i/ien we also have: 

(iv) ^ ^ 7 = 2|F | 2 |sym(g F - 1 )- 1 | 2 . 

(v) Z 2 ^ 7 = 2|F | 2 |sym(g F - 1 )- 1 | 4 . 

The proof of Lemma 12.31 will be given in section [5j We now remark that in the 
proof of the main Theorem 11.11 Lemma 12.11 will be used to the operators M = Mx = 
DqZ(F ,Qq), at finitely many spacial points X G Q, where F = W^o(X) an d Qo = 
V£i(X). It is clear that when the initial data £ 0; £1 with regularity (11.101) satisfy 
det V£ > (or the two conditions det V£ > and det sym(V£i(V£o)~ 1 ) 7^ whenever 
required) then the constants 7 in Lemma 12.31 have a common upper and lower bounds, 
independent of X. Therefore, Lemma [2.11 and the estimate ( I2.3P may be used with a 
uniform constant C Pi u, also independent of X. 

3. The main a-priori estimate 
Given £ 0) £1 as i n Theorem 11.11 let £1 G W 2,1 {VL x R + ) be the solution to: 

r (eOt-Ae 1 = o inoxt + , 

l & = on dVt x R + , 

I ^1(0, -) = ^1 in ft, 

Define the extension £ of £0, so that <9 f £ = £1: 

(3.1) e(t,x)=£ (x)+ f £i(s,x)ds. 

Jo 
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By continuity, it is clear that: inff2 X (o,T) det V£ > for T sufficiently small. We define: 

D = D(T) = ||^,V 2 6||L p{ nx(o,T)) 

and note that: 

(3.2) lim D(T) = 0. 

Lemma 3.1. Let £ ? £i be as ^ n Theorem \l.l\ and assume that: 

for every AJM holds with M = D Q Z(VZo(X),V£x(X)). 

Let £ G Wf } ' 2 {yt x (0,T )) with & G Wp 2 ' 1 ^ x (0,T )) be a solution to the problem / TO) . 
/ TO)) . ( TO)) . JTjp , and denote: 

e = 6(t) = - i) tt , v 2 (e - e)tlk(Qx ( o,T)), 

where £ is as in A3. Then, there exists T 00 < T and a constant C , both depending 
only on £o and £i (and, naturally, on Q andp), such that for every T < Too we have: 

(3.3) < C(T X I V + D + (T 1 ^ + D)G + 2 + 4 ). 
In particular: 

(3.4) 0(T) < C VT < T 00 . 

Before we give the proof of the lemma, we gather below some standard inequalities 
that will be frequently used for different functions: u defined on Q x (0,T), and w 
defined on Q. We always assume that T < 1. 

(3.5) ||Hlw|(fi) < C Pt n\\Aw\\ Lp (si) when w\ dn = 0, 



sup \\u(t, -)\\ 2-2 /Pm < C Pi n (\\ut- Au|| Lp( nx(o,T)) + \\u(0,-)\\ 2-y P(n) 

(3.6) *e(o,T) p \ > \ P \ > 

when n| 9 nx(o,T) = 0, 

(3.7) ||w|| ioo (n) < C P) q in fact: ||w||c°.«(n) < C a ^ p ^\\w\\ w i-i/ P ^ n y 



(3.8) ||Vw|| CQ ,«/2 < Ca,p,u\\u\\ w 2,\ Ux ( 0tT ))- 

The inequality ( 13. 5ft is the usual elliptic estimate [H], and ( 13. 6p is the parabolic estimate 
from [H] . The Morrey embedding gives ( 13.71) for p > n + 2 [H] , while (13.81) follows from 
the embedding VWp 2,1 (ft x (0,T)) C ^(U x (0,T)), also valid for p > n + 2 [17]. 
We stress that the constants C in all the above bounds are universal, i.e. they are 
independent of T. Additionally, the dependence of C in (13.81) on U C Q is uniform 
for any family of domains which are uniformly bilipschitz homeomorphic to each other 
after appropriate dilations. 



|V 2 u|| Lp (nx(0,T)) 




T 



Q JO 



V \ 1 /'P 

dt dX 



VISCOELASTICITY WITH PHYSICAL VISCOSITY 
We further remark the following simple bound: 

V 2 u t ds + V 2 u(0 

(3 9) / r r T \ 1/p 

< T VpM T p/p' J \V 2 u t \ p dtdX\ +T 1 / p ||V 2 n(0,-)||L p( n) 

= T 1 ^ (||W|U p( nx(0,T)) + ||V 2 M (0,-)||L p( n)) • 
Let now £ and £ be as in Lemma [3. 11 Using (13.61) to — we obtain: 

sup ||(£ -!)*(*, -)lk(n)+ sup ||V(£ -£)*(*, -)lk(n) 

(3.10) * e(0 ' T) <e(0,T) 

< C (II - Ott|U P (nx(o,T)) + l|V 2 (£ - Ot|U P (nx(o,T))) < C@, 

and consequently: 

(3.11) ||(£-£) t |k(nx(o,T))+ ||V(e-Otk(nx(o,r)) < CT^e. 
By ([3Z7]), dSSD used to £ - £, and (El), ( l3"TTTj) we get: 

||V(e-|)||L 00 (nx(o,T))= sup ||V(e-0(*,OIU=c(n) 

te(o,T) 

(3.12) < sup \\V(t-t)(t,-)\\ w ,- l/P 

te(o,T) p v ' 

< C(||(e-aiU p( nx(o,T)) + ||V 2 (e-0|U p( nx ( o,T))) < CT^e. 
Likewise, using (13. TP and (13.61) to — we directly obtain: 

(3.13) || (f - Otll^cnxCCT)) + || V(f - OtllL^axM) < CQ. 

In all the above inequalities ( 13 . 1 0[) - (13. 13j) . we write = O(T). The constant C 
depends only on the initial data of the problem £ , £i ( m addition to its dependence on 
O and p) . 

Proof of Lemma 13.11 

We will always assume that T < 1. Note that for T < Tq sufficiently small, the 
constraint (ll.4p is a consequence of the same constraint on the initial data £o, by 
continuity. Likewise: 

(3.14) \\DZ{Vl Vi t ),D 2 Z{Vl V{ t ),D 3 Z{Vl VZt)\\ Loo (tix(o,T)) < C. 
1. The system (11. ip can be rewritten as: 

(e - Ou - div (z(vf , V6) - z{yl v&)) = div (dw(vc)) + div (z(yl v&)) - e« 

and further, it has the form: 

(3.15) ^-0u-diy(D Q Z(VlVCt)V(C-0t) = F[^], 
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where: 



(3.16) 

F[£, £] = div (DW(VC)) + div (2(V£, V&)) - & 
+ div(Z(Vf,Vf t )-Z(Vf,V&)) 

+ div (z(ve, V6) - z{vl V6) - £> Q ^(ve, v&) v(e - 

= div (W(V0) + div (Z(Vf , V&)) - & 

+ div (L>^(ve,ve t )v(e-o) 

+ div (7 (1 - S )I^( S V£ + (1 - S )V£, V&)(V(£ - ® V(f - 0) dsj 



+ div ( y (1 ~ ^)^Q Q ^(Ve, + (1 - s)V&)(V(£ - 0* ® v(e - 00 
We shall now prove the bound: 

(3.17) \\F[£, e]|U P (nx(o,T) < C (T 1 ^ + D + (T 1 ^ + ,0)6 + 6 2 + 6 4 ) . 
By (J33D and (l3~T2|) it follows that: 

||div(W(VO)lk(nx ( o,T)) < \\D 2 W(V0\\l p (f2x(0,T)) 

<(\\d 2 w(vo\\l^ + c\mc-0\\l 



(3 ' 18) • (HV 2 (£ - 0t\\ Lp + \\V%\\ Lp + ||V 2 £o|lMn)) 

< C(l + T 1/P Q)T 1/P {1 + Q + D)< CT 1/P {1 + 6)(1 + 6 + D). 

Using (13.141) and (13.91) to £, we obtain: 



||div (Z(Vf,V&)) IU p( nx(o,T)) 
(3.19) < ||£>Z(V£ Ve t )|U oo( nx(o,T)) (II V 2 f || Lp(n x(o,r)) + ||V 2 &|| Mnx(0 ,T))) 

< C (T^|| V 2 6||l p + T l ' p \\ V%\\ Lim + \\V%\\ Lp ) < C(T^ + D). 
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By (EH, (EH, (EUD we get: 

(3.20) 

||div (D F Z{Vl V6) V(e - 0) IU P (nx(o,T)) 

< \\D 2 FF z{ylwi t )\\ Loo (||v 2 el M cx ( o,T)) + ||v 2 6||L p( nx(o,T))) l|v(e-OI| ioo 

+ || J D F Z(Ve,V6)||L oo ||V 2 (e-0l|L p( nx(0,T)) 

< c (i + ii v(e - &UO (ii v 2 elL P + ii v%\\ Lp + ii v 2 (^ - mw) II v(e - OIU. 

+c(i+iiv(e-e)tiio iiv 2 (e-oiu P 

< C(l + 0)(T 1/p + + D)T 1/p Q + C(l + 0)T 1/p 

< CT 1/p (l + + ^)0(1 + 0). 

and: 
(3.21) 

||div ( [ (1 - S )L> 2 ^( S V£ + (1 - s)V£, V&)(V(£ - ® V(£ - 0) ds) \\ Lp{ n xi o,T)) 



< sup ||div (D 2 FF Z(sV£ + (1 - s)V£, Ve*)(V(e - I) ® V(£ - 0)) lk(nx(o,T)) 

s€[0,l] 



< sup 

se[o,i] 



ii^( S ve+ a - S )ve, voikco (iiv 2 ^iil p + nv 2 (e-oiu P + nv 2 6iuj nv(e-oiiL 



+ || J D 2 z( s ve + (i- s )ve,V6)IU oo ||v 2 (e-0IU P ||v(e-0llL o 



< C (1 + || V(£ - 01k, + || V(e - OtllL.) (^ 1/P + © + D)T 2 / p Q 2 

+ C (1 + || V(£ - oik- + 1| v(e - T 2 / p 2 
<CT 1/p (l + + £)0 2 (l + 0). 

In the same manner, we see that: 
(3.22) 

||div (j (1 - s)D 2 QQ Z(Vl aVit + (1 - s)Vit){V{i - i) t ® V(£ - |)t) IU p( nx(o,T)) 

< C (1 + ||V(£ - &||lJ (||V 2 e|| Lp + ||V 2 (£ - t \\ Lp + ||V 2 6IUJ ||V(£ - aiiL 

+ c (i + ii v(e - e) t |iz.) ii v 2 (e - e)*iujiv(e - Otiko. 

< C (1 + 9) (T 1 ^ + + D)0 2 + C (1 + 0) 

< C(T 1 / p + + D)0(l + 0) 2 . 
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Combining ( 13381) - ( l3~22l) . the bound (13~T71) follows if only T < 1, ensuring £>(T) < 1 
by dS2D. 



2. We will now work with the localizations of the system (I3.15p . Let {B k }^ =1 be a 
covering of Q by a finite number N = N(r) of balls B k = B(X k , r) with centers X k G Q 
and radius r < 1. This family of coverings (parametrized by r) should such that all the 
sets 2Bk n Q are uniformly bilipschitz homeomorphic to each other after appropriate 
dilations and that the covering numbers of {2B k D £l}k are independent of r. 

Let 7Tfc : M n — > [0, 1] be smooth cut-off functions satisfying: n k = 1 on S^, and n k = 
on M n \25 fc where 25 fe = £(X fc , 2r), and ||V a 7r fc || Loo < CH a L After multiplying fl3TT5l) 
by 7Tfc, we obtain: 
(3.23) 

(vr fe (e - |))„ - div(D Q Z(VUXk), Va(^))V(7r fc (^ - |)0) = 7T fc F[£,£] + 
where: 



- div (D g 2:(VeoM, V6M)((£ - e)t ® Vvr,)) . 
We shall now prove the bound: 

||G^,£]|| ip(2SfcX (o,T)) <C(^ + T a/2 )||7T fc V 2 (e " Ot\\L p (2B k ,(0,T)) 

( 3 - 24 ) „, 1 



c(i + -)(T 1 / p + D)e(i + e). 

Using (I3.1ip . we obtain: 

|| (D Q Z(VZo{Xk),VZi(Xk)MZ - Dt) V7r fe || M2BfcX(0 , r)) 
(3-25) (7 C / 

<^l|v(e-0*lk<^e 2 . 

Likewise: 

||div (Dq^CV&W, VCl(X fc ))((e - Ot ® V7T fc )) ||L p ( Mh x(0,D) 

(3-26) n - C C , 

-\Mt-^ Lp + -M-0th P <- 2 T^e. 
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Finally, by 1ETI3]> . (13T9D . and (EHD we have: 

(3.27) 

||7r fc div([D Q z(vt, V6) - Aj^vCoM, v&(x fc ))]v(£ - £)*) H^b^d) 

< c(IIv 2 £1l p + ||v 2 6IU P ) ||v(£-£) t |k. 

+ hajZ(v£vo -^(veora,vei(x fc ))ii^ (2BfcX(0ir)) ii7r fe v 2 (e-e)*iu P ( 2 B fe x ( ( 

< C(T 1/p + .0)0 

+ c (nve - vto(x k )\\ Looi2 B k x(o,T)) + live* - vei(x fc )iu oo( 2B fc x(o,T) ) ) • 

• IkfeV '(£ - 0t|U p (2S fe x(0,T)) 

< C(T^ + D)0 + C(r Q + 2^)11^11^.1 Ik^Ce - £) t |U P (2B fe x(o,T)). 

Combining Q3.26P - Q3.27P and noting that ||^|| w -|.i(n X (o,T)) ^ ^(ll£o||wa(n)+||6|| w £_ a /p (n) ) : 
for T small enough we conclude (13.241) in view of (13.21) . 

We now use Lemma 12.11 to the problem (13.231) . i.e. we set ( = 7Tk(C, — £)t, M = 
DQZ(Vl;o{Xk),Vl;i(Xk)) where U = 2B k D Q is the uniform constant from the as- 
sumption (II. lip . Indeed, it is easy to notice that if (12.11) holds for some set U then it 
holds with the same constant 7 on every open subset U\ C U. By (12.31) we now obtain: 

IKfc(f - Ott, V 2 (7T fc (e - Ot)||i p (2B fe x(0,T)) < C r |kfci r [^|],G fc [^,^]|| £l ,(2B fc x(0 I T))- 

Summing over finitely many k : 1 . . . N, we get in view of (13.241) : 

II (e " 0«,V 2 (£ - OtlU p( nx(o,T)) < ^(r a + T Q / 2 )||V 2 (£ - $ t \\ Lp (nx { o,T)) 

+ cn i '*((i + 1/^(^ + 0)6(1 + 0) + iiF[e,e]ii M nx(o,T))), 

where, again, C depends only on the covering number of {B k } k=1 , on £, p and Q, but 
not on r, N, T or 0. Consequently, for r and T sufficiently small, we arrive at: 

< CN 1/p (l + ^ (T 1/p + D + (T 1/p + D)0 + 2 + 4 ) 

in virtue of fj3.17j) . This concludes the proof of (13. 3p . 

3. To prove ( 13.41) . consider the functions: 

g{Q) = Q and # e (0) = C(e + e© + 2 + 4 ), 

where C is a given constant and e > is a small parameter. 
Clearly, g(0) < g e (0) for every e. Take now: 

< 3 28) £<min <T^'4Z7' 1} 
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and let 9 e (4Ce, ^) with 6 < 1. Then: max{Ce, Ce 2 9 , CG 2 , C9^} < ^ and 
hence g(Q ) > g e (0o)- 

Taking now Too so small that, in addition to other requirements imposed in the course 
of the proof, e = T x ' v + D satisfies (I3.28p . we obtain that for every T e [0, T 00 ) the 
quantity 6(T) must stay below 9 , in virtue of continuity of the function T H- 6(T) 
and 6(0) = 0. This ends the proof of H3.4[) and of Lemma [37171 ■ 

4. A proof of Theorem 11.11 

We only outline the proof of Theorem ll.il which is standard, and we point to its most 
important steps. Let £ be as in (13.11) . Recall that the system (II .ip can be rewritten as: 

(4-1) (Z-Ott-div(D Q Z{VtM t )V(t-£) t ) =F[£,|], 

where the right hand side £] is given in (I3.16p . We shall seek a solution £ as the 
fixed point of the operator: 

T(£ — £) = £ — £5 £isa solution to: 

(4-2) (C - On - div (D Q Z{Vl V&) V(£ - |) t ) = F[£, £], 

in the Banach space: 

E n , T ={ue L P (Q x (0, T)); u(0, •) = 0, u t (0, •) = 0, «| M x(o,t) = 0, 
(4-3) 1 

%,V\eL p (Ox(0,T))j, 

equipped with the norm: 

IMUn.T = ©MCO = \\ u tt, V 2 n t || LP(nx (o,T). 

1. Following calculations as in the proof of Lemma [3. 1\ it results that: 

Vf - £ e E U>T F[l £] e L p (fi x (0, T). 

2. Integrating (14. 2 p against (£ — £)t on f2 x (0,T) and using the estimate (14. 5 p in 
Lemma [4.11 below with ( = (£ — £) t , we obtain: 

sup ii(£-o*(v)IIL ( n) + llv(e-o t |ii 2(nx(0) T)) 

<e(o,r) 

< C T / dx dt + C||(£-O t ||l 2(nx(0iT)) 

Jo jq 

<C/ T /" F%Z}dxdt + CT sup ||(£-&(V)IIL(0)> 
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which implies the following energy estimate, for T small: 

(4.4) sup ||(£-£)t(V)ll L(nx(o,T)) 

te(o,T) Jo Jn 

In virtue of (I4.4p . the Galerkin construction of the approximants: 

N 

(6v-£)t = a&(t) wifc), 

k=l 

where {wi}^ is an orthonormal base of W^i^i), yields existence of a weak solution 
£ - C = lim^oo^Tv - of the problem g^D, with: (f - |) f £ ^((0, T), L 2 (Q)) and 
V(e-0t e£ 2 (fix (0,T)). 

3. A modification of arguments in section [3] implies that the weak solution £ is 
actually regular in the class determined by (14. 3p . i.e: 

Moreover, for every small e > 0: 

if G[i- £](T) < e then 6[£ - £](T) < e. 

4. In now suffices to show that the map T is a contraction in some ball B e C Eq^- 
This is done by applying methods of (E]) to the system: 

(f i - 6)« - div (D Q z(vt, vcd v(6 - 6)0 = ^[£1, e] - Ffe, e], 

where T(£i — £) = 6 — £. For e > sufficiently small it follows that: 

e[6-e 2 ](r)<^©[fi-6](r), 

which completes the proof. ■ 

The key role above was played by the following estimate: 

Lemma 4.1. Let T < T be sufficiently small and assume that Z satisfies U.ll\) . Then 
for every ( £ (0 x (0, T)) snc/i £/ia£ £(0, •) = and C|<9Qx(o,t) = 0, there holds: 

( 4 - 5 ) II VCllL ( nx(o,T)) < 47 T / A>2(V£ V6)VC : VC d* dt + C||C||i 2(f , x(0 . r)) , 
rt/i constant C independent of £. 

Proof. Consider a covering {-B^}^ of by a finite number iV = iV(r) of balls = 
B(Xk, r) with centers Xf. £ Q and radius r > 0. This family of coverings (parametrized 
by r) should be such that their covering numbers are uniform in r. Let {TTk}k=i be a 
partition of unity subject to {B k }. 
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For a fixed t e (0, T), with a slight abuse of notation, we shall still write ( = £(£, •) G 
W^iVt). By ffTTTTj) it follows that: 

f (D Q Z(V^^i)VC:VC) 
Jn 

N 

(4.6) =£/ (^(veo™,veira)v(7r; /2 c):v(7r; /2 o)dx + ^ / E fe [e,el 
>-Eii v ^ /2 oiiW) + E/ 



k=l 

where we accumulated the error terms in: 

1/2,-,* . 1/2 



- <A?Z(V£ , V6)V(7r;/ 2 C) : V(7r;/ 2 C)) 

+ ([D Q Z(VCo, V6) " A^V^X*), V{ 1 (X fc ))]V(7r; /2 C) : V^C)) 

Hence: 

| /" £ fc [£,£] dx\< \ [ (D Q Z(V^^iH /2 ^C ■ (C® Vvr; /2 ))| 

+ | / (D Q Z(Veo,V6)(C®V7rf ): V(7rfO)l + Cr||V( 7 ri /2 C)lli 2 ( Bfe ) 



< 



Cr||Cll^(B fc )llClU 2 (B fc )+C , r||V(7r; /2 C)ll! 2( i ?fc) 



where C is a universal constant depending only on the initial data and Z, while the 
constant C r depends on the covering {B r }. Taking r small, so that Cr < 1/(27), by 
(14. 6p we now arrive at: 



f (A^v^veovc : vc> 



N 

> 

2^ 

k=l 

N 



1 

2^En v (^ /2 c)HW)-aiicii^ ( o)iiciu 2 (n) 



>^Eii^ /2v ciiL( Bfe )-aiicii^(o)iiciu 2 (o) 
~ fc=i 

> ^l|VCIlL(n) - ^ (ellVfllLpi) + ^IIVCIIL(O) i , 
by Young's inequality. With e sufficiently small, it yields: 

(4-7) iivciiLm ^ / <A?2(V£ , veovc : vc> + C||C|| 



I ia(«)- 
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Integrating in t, we eventually arrive at: 



(f2x(0,T)) 



l|VC|lL(nx(o,D) < 37 f f (^(VCo, Vei)VC : VC> dx dt + C\\C\\% 
Jo Jn 

< 3 7 f [ (D Q Z(Vl V6)VC : VC) dx dt 
Jo Jn 

+ CT||VC|| L 2 (nx(0,T)) + C||CII L 2 (nx(0,T))) 
which for T small enough implies (14.51) . 



5. A proof of Lemma 12.31 

1. To prove (i), note that DqZq(F ,Q )Q = 2F sym(F T Q) so that: 

VQ e R nxn (D Q Z' \F , Q )Q : Q) = 2(sym(F T Q ) : F*Q) = |sym(F T Q)| 2 . 
Take ( G ^(f^IR™) with trace on the boundary d£l. We have: 

f |VC| 2 < |F "^| 2 f |V(F r C)| 2 <2|F ^| 2 / |symV(ifC)| 2 
Jn Jn Jn 

= IF- 1 ' 7 ] 2 [ (D Q Z'J(F ,Q )V( : VC), 
Jn 

where we applied Korn's inequality to the map x i— > Fq£(x). 

2. To prove (ii), observe that DqZq(F ,Q )Q = 2(det F )sym(QF Q ~ 1 )F ~ 1,T so that: 

(D Q Z' (F ,Qo)Q : Q) = 2(detF )|sym(QF - 1 )| 2 . 
Then, for any test function C as above, we have: 

f |VC| 2 <|F | 2 / |(VC)F " 1 | 2 dx=|F | 2 f |V(Co(F " 1 y))| 2 (detF " 1 )dy 
Jn Jn Jf u 

< 2|F | 2 (detF - 1 ) f |sym((VC)oF ( 7 1 )F c r 1 )| 2 dy 

(5.1) Jp ° n 

= 2|F | 2 (detF - 1 )(detF ) / \sym((VC)F^)\ 2 dx 

Jn 

= |F | 2 (detF )^ / (D Q Z' (F ,Q )V{ : VC), 
where we applied Korn's inequality to the map y h-> ((Fq 1 i/) on the open domain F Q. 
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3. To prove (iii) - (v), observe that: 



2 m 



[D Q Z m (F ,Q )Q :Q) = ( ^(sym(Q ^ 1 ) J sym(QF - 1 )(sym(g F c r 1 ) 2m ^ : QF^ 1 

3=0 
2m 

A j BA 2m - j : QFq 1 



3=0 

where we denoted: 

A = symiQoF,- 1 ), B = sym{QF^). 
Since the matrix Y^j2oA^BA 2m ~-' is symmetric, it follows that: 

2m 

(D Q Z m (F ,Q )Q : Q) = (Y, A3BA2m ~ 3 ■ B ] 

j=0 

Let ( be a test function as in Lemma 12 . 11 By calculations similar to (15. ip we get 

[ |VC| 2 < ^|F | 2 f (D Q Z (F ,Qo)VC : VC), 
Jn A Jn 

proving (iii). To prove (iv), we compute: 

(DqZ^Fq, Q )Q : Q) = (A 2 B : B) + (ABA : B) + (BA 2 : B) 

= (AB : AB)(BA : AB) + \AB\ 2 = 2(sym(AB) : AB) + \AB\ 2 

= 2|sym(AB)| 2 + \AB\ 2 > \AB\ 2 . 
Therefore, by calculations similar to ( 15. ip : 

f |VC| 2 <2|F | 2 f |sym((VC)F - 1 )| 2 <2|F | 2 |A" 1 | 2 / \AB\ 2 
Jn Jn Jn 

< 2|F | 2 |sym(g F " 1 )" 1 | 2 f (D Q Z 1 (F ,Q )V( : VC). 

Jn 

Finally, in order to prove (v) we derive: 

(D Q Z 2 (F , Q )Q : Q) = (A 4 B + A 3 BA + A 2 BA 2 + ABA 3 + BA 4 : B) 

= \A 2 B + ABA\ 2 + \A 2 B\ 2 > \A 2 B\ 2 , 
which, in the same manner as in (15 .2p yields: 

2 <2|F | 2 |yT 2 | 2 f \A 2 B\ 2 



(5.2 



[ IVCP 

Jn 

<2|F | 2 |sym(g i 71 o" 1 )" 1 | 4 / (D Q Z 2 (F ,Q )V( :V(). 

Jn 



The proof of Lemma [2.31 is done. 
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